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Abstract 

Programs  for  integrated  circuit  layout  typically  have  two 
phases:  placement  and  routing.  The  router  should  produce  as 
efficient  a  layout  as  possible,  but  of  course  the  quality  of  the 
routing  depends  heavily  on  the  quality  of  the  placement.  On 
the  other  hand,  the  placement  procedure  ideally  should  know 
the  quality  of  a  routing  before  it  routes  the  wires.  In  this  talk 
we  present  an  optimal  solution  for  a  practical,  common  version 
of  this  placement  and  routing  problem. 

River  routing  is  the  problem  of  connecting  in  order  a  set  of 
terminals  a\ , . . .  ,an  on  a  line  to  another  set  b\ , . . . ,  bn  across  a 
rectangular  channel.  Since  the  terminals  are  located  on  modules, 
the  modules  must  be  placed  relative  to  one  another  before  rout¬ 
ing.  This  placement  problem  arises  frequently  in  design  sys¬ 
tems  like  bristle- blocks  where  stretch  lines  through  a  module  can 
effectively  break  it  into  several  chunks,  each  of  which  must  be 
placed  separately.  In  this  talk  we  shall  present  concise  neces¬ 
sary  and  sufficient  conditions  for  wirability  which  are  applied 
to  reduce  the  optimal  placement  problem  to  the  graph-theoretic 
single-source-longest-paths  problem.  By  exploiting  the  special 
structure  of  graphs  that  arise  from  the  placement  problem  for 
rectilinear  win  \g,  an  optimal  solution  may  be  determined  in 
linear  time. 

Key  words  and  phrases:  analysis  of  algorithms,  graph  theory,  longest  paths,  placement  and 

routing,  river  routing,  VLSI  layout. 
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1 .  Introduction 


River  routing  is  a  special  routing  problem  which  arises  often  in  the  design  of  integrated 
circuits,  and  it  has  been  shown  to  be  optimally  solvable  in  polynomial-time  for  many  wiring 
models  (see  in  particular  [Tompa)  and  [Dolev  et  a/.]).  In  this  paper  we  demonstrate  that  the 
placement  problem  for  river  routing  is  also  polynomial-time  solvable. 

The  general  character  of  the  placement  problem  for  river  routing  is  illustrated  in  Figure  I. 
Two  sets  of  terminals  a\% . . . ,  an  and  b\9 . . . ,  bn  are  to  be  connected  by  wires  across  a  rectangular 
channel  so  that  wire  i  is  routed  from  at  to  b*.  The  terminals  on  each  side  of  the  channel  are 
grouped  into  chunks  which  must  be  placed  as  a  unit.  The  quality  of  a  legal  placement-  one  for 
which  the  channel  can  be  routed — can  be  measured  in  terms  of  the  dimensions  of  the  channel. 
The  separation  is  the  vertical  distance  between  the  two  lines  of  terminals,  and  the  spread  is  the 
horizontal  dimension  of  the  channel. 
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Figure  1:  Two  sets  of  chunks  on  either  side  of  a  rectangular  channel. 

Terminal  at  must  be  connected  to  b%  for  i  =  1, . . . ,  10. 

The  wiring  model  gives  the  constraints  that  the  routing  must  satisfy.  Although  our  results 
can  be  generalized  to  include  a  variety  of  wiring  models  (see  Section  5),  we  concentrate  on  the 
(one 'layer)  square' grid  model.  Crossovers  are  disallowed  in  the  square-grid  model,  and  all  wires 
must  take  disjoint  paths  through  the  grid. 

The  placement  problem  for  river  routing  arises  often  during  ordinary  integrated  circuit  design. 
A  common  instance  is  when  the  terminals  of  one  or  more  modules  are  to  be  connected  to  drivers. 
The  various  independent  “chunks”  are  the  modules,  which  lie  on  one  side  of  the  channel,  and  the 
drivers,  which  lie  on  the  other. 

A  more  interesting  manifestation  of  the  placement  problem  occurs  in  the  context  of  such 
design  systems  as  bristle-blocks  fJohannsen)  and  DPL/Daedalus  (Batali  et  a/.j.  These  systems 
encourage  a  designer  to  build  plug-together  modules  so  that  the  difficulties  associated  with  general 
routing  can  be  avoided.  A  designer  may  specify  stretch  hues  which  run  through  a  module  and 
allow  the  module  to  be  expanded  perpendicular  to  the  stretch  line,  as  demonstrated  in  Figure 
2.  When  two  independently  designed  modules  are  plugged  together,  stretch  lines  permit  the 
terminals  to  be  pitch  aligned ,  that  is,  the  distances  between  pairs  of  adjacent  terminals  are  made 
to  match  the  distances  between  their  mates,  and  routing  is  avoided  because  the  Reparation  of 
the  channel  is  zero.  Unfortunately,  this  approach  may  not  succeed  unless  stretch  lines  are  put 
between  every  pair  of  adjacent  terminals.  The  stretch  lines  may  not  only  disrupt  the  internal 
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structure  of  the  modules,  but  the  cons*  quence  may  be  an  inordinate  amount  of  stretching  that 
leaves  the  channel  with  a  large  spread. 


The  other  extreme  is  to  forego  stretching  altogether  and  river  route  between  the  terminals. 
But  the  cost  may  still  be  large  if  a  large  separation  is  required  in  order  to  achieve  a  routing. 
A  reasonable  compromise  is  to  place  stretch  lines  where  it  is  convenient,  and  then  do  a  little 
stretching  and  a  little  routing.  Determining  how  much  of  each  to  do  is  exactly  the  optimal 
placement  problem  for  river  routing. 


Figure  2:  A  module  before  and  after  stretching  (courtesy  of  John  Batali). 

The  remainder  of  this  paper  demonstrates  that  optimal  solutions  to  the  placement  problem 
ran  be  achieved  efficiently.  Section  2  gives  a  concise  necessary  and  sufficient  condition  for  a 
channel  to  be  routable  in  the  square-grid  model.  Section  3  shows  that  the  form  of  this  condition 
allows  the  placement  problem  to  be  reduced  to  the  graph- theoretic  problem  of  finding  the  longest 
paths  from  a  source  vertex  to  all  other  vertices  in  a  graph.  Based  on  this  problem  reduction, 
a  linear-time  algorithm  for  optim  1  placement  is  given  in  Section  4,  Section  5  shows  that  the 
algorithm  extends  to  wiring  models  other  than  the  square  grid  model,  but  its  performance  depends 
on  the  particular  wirability  conditions  for  the  model.  Section  6  discusses  the  application  of  our 
results  to  other  routing  situations  and  suggests  further  placement  problems. 
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2.  Necessary  and  Sufficient  Conditions  for  Wirability 

To  demonstrate  the  results  of  this  paper,  we  adopt  an  extremely  simple  wiring  model:  the 
(one-layer)  square-grid  model.  All  wires  are  constrained  to  run  on  an  underlying  grid  of  integer 
lattice  points.  The  terminals  a j , . . . ,  an  and  bi , . . . ,  bn  occupy  grid  points  on  opposite  sides  of  the 
channel.  No  two  wires  may  occupy  a  single  grid  point  which  enforces  unit  separation  of  wires. 
Figure  3  shows  a  solution  to  the  problem  of  Figure  1  using  this  model. 


Figure  3:  A  possible  solution  to  the  problem  in  Figure  1  for  which  the 
separation  is  5  and  the  spread  is  27. 

In  order  to  establish  constraints  on  wirability  in  this  model,  consider  a  straight  line  segment 
drawn  from  to,  but  not  including,  (12,2/2)-  We  ask  the  question,  “How  many  wires  can 

cross  this  line?”  A  simple  analysis  shows  the  answer  is  max(|i2  —  X\ |,  I2/2  —  y i  ()•  Without  loss  of 
generality,  assume  the  situation  is  as  in  Figure  4,  and  look  at  the  grid  points  immediately  below 
the  line,  that  is,  { (xf  —  ^i)J)  |  x  =  x\t . . . ,  Zg  —  1  }.  Any  wire  crossing  the  line 

must  perforce  occupy  one  of  these  grid  points,  and  therefore  the  number  of  such  wires  is  bounded 
by  the  cardinality  of  this  set. 


(*2,1/2) 


Figure  4:  The  number  of  wires  crossing  the  half-open  line  segment  is  at 
most  the  number  of  grid  points  immediately  below  the  line. 

Let  us  now  turn  to  the  river  routing  problem  and  examine  how  this  constraint  can  be  brought 
to  bear.  Let  ai,...,an  denote  both  the  names  of  the  terminals  at  the  top  of  the  channel  and 
their  z-coordinates,  and  let  the  same  convention  hold  for  the  terminals  bi,...,b„  at  the  bottom 
of  the  channel.  Consider  a  half-open  line  segment  drawn  from  terminal  a?  to  terminal  bx  as  shown 
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ir  Figure  5.  The  j  —  i  wires  emanating  from  a,,...  ,  aj_i  must  all  cross  this  line,  and  similarly 
f  ^r  a  line  drawn  from  b3  to  ax.  In  order  for  a  channel  with  separation  t  to  be  routable,  therefore, 
ic  must  be  the  case  that 

max(a,  —  b%f  t)  >  j  —  i  and  max(d;  —  att  t)  >  j  —  i  (1) 

for  1  <  i  <  j  <  n. 


j  —  i  wires 


✓ 

✓ 

/ 

/ 

s 

* 


bx  bx+1  ...  6,_i  bj 


Figure  5:  The  j  —  i  wires  from  a%, . . . ,  must  cross  the  dashed  line 
between  b%  and  ay 

Although  Condition  (1)  is  a  new  condition  for  wirability,  the  analysis  that  leads  to  it  is 
essentially  the  same  as  that  in  [Dolev  et  al. ]  and  represents  previous  work  in  the  field.  A  more 
compact  condition  exists,  however,  which  is  equivalent: 

at+t  —  bx  >  t  and  bx+t  —  a%  >  t  (2) 

for  1  <  i  <  n  —  t.  The  channel  is  always  routable  if  t  >  n. 

Condition  (1)  implies  Condition  (2)  because  Condition  (2)  is  a  refinement  of  Condition  (1). 
For  the  opposite  direction,  suppose  first  that  j  —  i  <  t;  then  max(a,  —  6*,  t)  >  t  >  j  —  i.  If 
j  —  i  >  t<  on  the  other  hand,  then 

aj  ~~  h  =r  o  ~~  h 

>  *%+t  —  +  —  i  — 

>  t  +  (j  —  i  —  t) 

since  \  >  a*  +  1  for  all  1  <  k  <  n.  Thus  the  two  conditions  are  indeed  equivalent. 

Figure  6  shows  a  simple  geometric  interpretation  of  Condition  (2).  The  condition  a%+t  —  bt  > 
t  means  that  a  line  with  unit  slope  going  up  and  to  the  right  from  b%  must  intersect  the  top  of 
the  channel  at  or  to  the  left  of  terminal  And  if  the  condition  fails,  terminal  b7  must  be  to 

I  he  right  of  a3  for  i  <  j  <  i  -f- 1  —  1,  that  is,  each  wire  from  a  a3  goes  down  and  to  the  right, 
which  can  be  shown  to  follow  from  the  fact  that  a3  +  \  >  a3  *f-  1.  (For  6t+t  —  a,  >  t  the  line 
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Figure  6:  Geometric  interpretation  of  ax+t  >  K  +  t  and  6v.^t  >  at  +  ** 

with  slope  — 1  going  down  and  to  the  right  from  a,  must  intersect  the  bottom  of  the  channel  at 
or  to  the  left  of  terminal 

This  geometric  interpretation  can  be  used  to  show  that  Condition  (2)  is  not  only  a  necessary 
condition  for  routabiiity  of  the  channel,  but  a  sufficient  condition  as  well.  In  fact,  a  simple  greedy 
algorithm  will  successfully  route  a  routable  channel.  Processing  terminals  left  to  right,  the  greedy 
algorithm  routes  each  wire  across  the  channel  until  it  hits  a  previously  routed  wire;  then  it  follows 
the  contour  of  the  opposite  side  until  it  reaches  its  destination. 

To  see  that  this  algorithm  works  given  Condition  (2),  wc  must  he  more  prerise  about  what 
paths  are  taken  by  the  wires.  Consider  without  loss  of  generality  a  block  of  consecutive  wires 
that  go  down  and  to  the  right,  that  is,  a,  <  ht  for  all  wires  in  the  block.  For  any  horizontal 
position  x  such  that  at  —  t  <  x  <  bx>  define 


r)t(x)  =  max(a,  —  x,  max  r). 

b,_r>  I 

The  path  of  wire  i  is  then  described  by  the  locus  of  points  (x  +  rjt(x),  J?»(x))  for  at  —  t  <  x  <  bx. 

A  geometric  interpretation  of  this  formulation  uses  the  same  intuition  as  was  given  in  Figure 
6.  The  line  with  unit  slope  drawn  from  (x,  0)  where  x  is  in  the  range  a%  —  t  <  x  <  bt  must  cross 
wire  i.  The  value  rjt(x)  gives  the  y-coordinate  of  wire  i  where  it  crosses  this  line  of  unit  slope. 
The  two-part  maximum  in  the  definition  of  Vx{x)  corresponds  to  whether  the  wire  is  being  routed 
straight  across  the  channel  or  whether  it  is  following  the  contour  of  the  bottom.  The  value  of 
rjx(z)  for  the  latter  situation  is  the  number  of  wires  to  the  left  of  wire  i  which  must  cross  the  line 
of  unit  slope. 

We  must  now  show  that  the  locus  of  points  for  a  wire  is  a  path,  that  the  paths  are  disjoint, 
and  that  they  never  leave  the  channel.  That  the  locus  of  points  is  indeed  a  path  can  be  seen 
by  observing  that  as  x  ranges  from  ax  —  t  to  bt,  the  initial  point  is  (o,,f  —  1),  the  final  point 
is  (6^0),  and  with  a  change  of  one  in  x  the  coordinates  of  the  path  change  by  a  single  grid 
unit  in  exactly  one  of  the  two  dimensions.  To  show  that  the  paths  are  disjoint,  consider  two 
adjacent  wires  i  and  i  +  h  and  observe  for  +  j  —  t  <  x  <  that  a,  —  x  <  *  t  —  x  and 
max&,„r>z  r  <  max*t4_,  _r>*  r,  and  therefore  rjx{x)  <  r?t+l(x). 

To  show  a  path  of  a  wire  never  leaves  the  channel,  we  demonstrate  that  rj»(x)  <  t  for  all  i 
and  x  in  the  associated  range.  It  is  for  this  part  of  the  proof  that  we  need  the  assumption  that 
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Condition  (2)  holds.  If  for  a  wire  i,  the  two-part  maximum  in  the  definitior  ol  rjx(z)  is  achieved 
">y  at  —  x,  then  ^(x)  must  be  less  than  t  because  x  >  a,  —  t.  Suppose  then,  that  the  two-part 
maximum  is  achieved  by  the  maximal  r  such  that  6t_r  >  x.  To  show  that  r  <  t,  w <;  assume  the 
contrary  and  obtain  a  contradiction.  But  since  >  bx—r  >  x  >  ax  —  t,  the  contradiction  is 
immediate  because  at  —  bx—t  >  t  from  Condition  (2). 


3.  The  Structure  of  the  Placement  Problem 

The  objective  of  a  placement  algorithm  is  to  set  up  a  routing  problem  that  is  solvable  and 
minimizes  some  cost  function.  Many  criteria  can  be  adopted  to  measure  the  cost  of  a  placement 
for  river  routing,  whether  in  terms  of  area  (total  or  channel)  or  some  other  function  of  spread  and 
separation.  A  plot  of  minimal  spread  versus  given  separation  reveals  that  the  region  of  feasible 
placements  may  not  be  convex  although  the  curve  is  guaranteed  to  be  monotonically  decreasing. 
(Figure  7  shows  the  plot  for  the  problem  of  Figure  1.)  Any  measure  of  placement  cost  that  is  a 
function  of  spread  and  separation  and  which  is  monotonically  increasing  in  each  of  spread  and 
separation  will  therefore  find  a  minimum  on  this  curve. 
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Figure  7:  The  curve  of  minimum  spread  versus  separation  for  the  example 
of  Figure  1. 

Thus  wo  content  ourselves  with  producing  points  on  this  curve,  that  is,  determining  a 
placement  which  achieves  the  minirnum  spread  for  a  given  separation  t ,  if  indeed  the  channel  is 
routable  in  t  tracks.  If  minimum  separation  is  the  goal,  for  example,  binary  search  can  determine 
f hr*  optimum  t  in  0(\e,t)  steps.  Since  the  algorithm  presented  in  the  next  section  determines  a 
placement  for  fixed  t  in  0(?i)  time  where  n  is  the  number  of  terminals,  and  since  the  separation 
need  never  be  more  than  n,  a  minimum-separation  placement  can  be  achieved  in  O(nlgn)  time. 
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F  >T  more  general  objective  functions  such  as  area,  the  optimum  value  can  be  determined  in  0(n2) 
t*me. 

We  now  examine  the  character  of  the  placement  problem  for  river  routing  when  the  separation 
t  is  given.  The  n  terminals  are  located  on  m  chunks  which  are  partitioned  into  two  sets  that 
form  the  top  and  bottom  of  the  channel.  For  convenience,  we  shall  number  the  chunks  from  one 
to  k  on  the  top,  and  k  +  1  to  m  on  the  bottom.  The  order  of  chunks  on  each  side  of  channel 
is  fixed,  but  they  may  be  moved  sideways  so  long  as  they  do  not  overlap.  For  each  chunk  i,  a 
variable  vt  represents  the  horizontal  position  of  its  left  edge.  Any  placement  can  therefore  be 
specified  by  an  assignment  of  values  to  these  variables.  We  also  add  two  variables  Vq  and 
to  the  sot  of  variables,  which  represent  the  left  and  right  boundaries  of  the  channel.  The  spread 
is  thus  Vm- |-i  —  Vq,  Figure  8(a)  shows  the  eight  variables  for  the  example  from  Figure  1. 
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(a)  Assignment  of  variables  to  chunks  and  channel  boundaries. 


Figure  8:  Representing  the  placement  constraints  as  a  graph  for  the  ex¬ 
ample  of  Figure  1. 

Since  the  relative  positions  of  terminals  within  a  chunk  is  fixed,  the  wirability  constraints  of 
Condition  (2)  can  be  reexprrssed  in  terms  of  the  chunks  themselves  to  give  placement  constraints 
that  any  assignment  of  values  to  the  v%  must  satisfy.  If  terminal  a,+t  lies  on  chunk  h ,  and 
terminal  lies  on  chunk  7,  the  constraint  a,+t  —  bx>t  can  be  rewritten  as  v ^  —  v3  >  r^,  where 
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rkj  reflects  t  and  the  offsets  of  the  terminals  from  the  left  edge  of  their  res  active  chunks.  The 
constraint  between  two  chunks  determined  in  this  way  will  be  the  maximal  constraint  induced 
by  pairs  of  terminals. 

Additional  constraints  arise  from  the  relative  positions  of  chunks  on  either  side  of  the  channel. 
For  each  pair  of  adjacent  chunks  i  and  i  +  1,  the  constraint  vt+\  —  v%  >  wt  must  be  added  to 
the  set  of  placement  constraints,  where  wt  is  the  width  of  chunk  t.  Four  more  constraints  are 
needed  which  involve  the  boundary  variables  and  um+i.  For  chunks  1  and  k  -F  1  which  are 
leftmost  on  the  top  and  bottom,  the  constraints  v\  —  vq  >0  and  Vk-\-i  —  vo  >  0  enforce  that 
these  chunks  lie  to  the  right  of  the  left  boundary  of  the  channel.  For  chunks  k  and  m  which 
are  rightmost  on  the  top  and  bottom,  the  relations  vm_j_ i  —  Vk  >  and  vm+\  —  Vm  > 
constrain  them  to  lie  to  the  left  of  the  right  boundary,  where  Wk  and  wm  are  the  widths  of  the 
chunks. 

Figure  8(b)  shows  a  placement  graph  which  represents  the  constraints  between  chunks  for  the 
placement  problem  of  Figure  1  where  the  separation  is  3  tracks.  A  directed  edge  with  weight  Ski 
goes  from  Vk  to  Vi  if  there  is  a  constraint  of  the  form  Vi  —  Vk  >  6kl-  For  example,  the  weight  of  1 
on  the  cross  edge  going  from  to  v2  is  the  maximal  constraint  of  ag  —  bG  >  3  and  ax o  —  67  >  3 
which  yield  u2  —  >  — 1  and  u2  —  U5  >  1  since  ag  =  t>2  +  5,  a  10  =  v2  +  6,  &6  =  vs  ~F 

and  67  =  -h  4.  The  side  edge  from  v+  to  V5  arises  from  the  constraint  that  chunk  4,  which  is 
5  units  long,  must  not  overlap  chunk  5. 

The  goal  of  the  placement  problem  is  to  find  an  assignment  of  values  to  the  vt  which  minimizes 
the  spread  y  —  Uq  subject  to  the  set  of  constraints.  This  formulation  is  an  instance  of  linear 
programming  where  both  the  constraints  and  the  objective  function  involve  only  differences  of 
variables.  Not  surprisingly,  this  problem  can  be  solved  more  efficiently  than  by  using  general 
finear  programming  techniques.  In  fact,  it  reduces  to  a  single-source-longest-paths  problem  in 
the  placement  graph.  The  length  of  a  longest  path  from  vq  to  corresponds  to  the  smallest 

spread  of  the  channel  that  complies  with  all  the  constraints.  The  placement  of  each  chunk  i 
relative  to  the  left  end  of  the  channel  is  the  longest  path  from  v0  to  vt.  If  the  placement  graph 
has  a  cycle  of  positive  weight,  then  no  placement  is  possible  for  the  given  separation. 

For  the  placement  problem  of  Figure  1  with  a  three-track  separation,  the  longest  path  from 
v0  to  v-i  in  the  placement  graph  (Figure  8)  is  Vo  —  vy  —  ^4  —  ^5  —  ^2  with  weight  13  which 
corresponds  to  the  positioning  of  chunk  2  in  the  optimal  placement  shown  in  Figure  9(a).  Figures 
9(b)  through  9(d)  show  optimal  solutions  to  the  placement  problem  of  Figure  1  for  separations 
t  =  4  through  t  =  6.  The  constraints  for  t  =  2  yield  a  cycle  of  positive  weight  in  the  placement 
graph,  and  thus  no  placement  is  possible  which  achieves  a  separation  of  only  two  tracks. 


4.  A  Linear-Time  Algorithm  for  the  Placement  Problem 

The  analysis  of  Section  3  showed  that  the  optimal  placement  problem  for  fixed- separation 
river  routing  was  reducible  to  the  single-source-longest-paths  problem  on  a  placement  graph. 
For  a  general  graph  G  =  (V',/?)  this  problem  can  be  solved  in  time  0(|V|  ■  |/r|)  by  a  Dellman- 
Ford  algorithm  (Lawler),  Better  performance  is  possible,  however,  due  to  the  special  structure  of 
placement  graphs.  This  section  reviews  the  Bellman-Ford  algorithm,  and  shows  how  it  can  be 
adapted  to  give  an  0(m)-time  algorithm  for  the  longest-paths  problem  on  a  placement  graph, 
where  rn  is  the  number  of  chunks.  Since  the  placement  constraints  can  be  generated  in  0(n) 
time,  where  n  is  the  number  of  terminal  pairs,  this  algorithm  leads  to  an  optimiil  linear-time 


8 


(a)  Separation  3,  spread  27. 


(b)  Separation  4,  spread  26. 


Figure  9:  Optimal  placements  and  routings  for  the  problem  of  Figure  1 
with  separations  ranging  from  t  =  3  to  t  =  6. 


a’gorithm  for  the  fixed-separation  placement  problem.  The  discovery  of  a  linear-time  algorithm 
represents  joint  research  with  James  B  Saxe  of  Carnegie-Mellon  University. 

The  linear-time  algorithm  is  a  refinement  of  the  standard  Bellman-Ford  algorithm  which  for 
each  vertex  vx  where  i  =  1, . . . ,  m  +  1,  iteratively  updates  the  length  X(vt)  of  a  tentative  longest 
path  from  v0  to  vt.  The  algorithm  initializes  X(u0)  to  zero,  and  all  other  X(u,)  to  —  oo;  then  it 
sequences  through  a  list  <f  of  edges,  and  for  each  edge  (vt)  v;)  with  weight  updates  X(vj)  by 

\(vj)  «-  max(X(vj),  6t,  +  X(v,)). 

The  list  £  of  edges  is  the  key  to  the  correctness  of  the  algorithm.  The  length  of  a  longest 
path  from  the  source  v0  to  a  vertex  v3  converges  to  the  correct  value  if  the  edges  of  the  path  form 
a  subsequence  of  the  list  £ .  (This  can  be  proved  by  adapting  the  analysis  of  [Yen].)  In  the 
normal  algorithm  for  a  general  graph  G  =  {V,E)}  the  list  £  is  \  V\  —  1  repetitions  of  an  arbitrary 
ordering  of  the  edges  in  E ,  which  ensures  that  every  vertex-diBjoint  path  in  G  beginning  with 
u0  is  a  subsequence  of  £ .  If  there  are  no  cycles  of  positive  weight  in  the  graph  Gy  then  from  t>o 
to  each  other  vertex  in  G,  there  is  a  longest  path  that  is  vertex-disjoint;  hence  the  algorithm 
is  guaranteed  to  succeed.  The  condition  of  positive- weight  cycles  can  be  tested  at  the  end  of 
the  algorithm  either  by  checking  whether  all  constraints  are  satisfied  or  by  simply  running  the 
algorithm  through  the  edges  in  E  one  additional  time  and  testing  whether  the  values  of  any  X(vt) 
change. 

The  list  if  is  also  the  key  to  the  performance  of  a  Bellman-Ford  algorithm.  For  the  general 
algorithm  on  an  arbitrary  graph  G  =  (F,  E),  the  length  of  the  list  is  (|F |  —  1)  •  |Z5|,  and  thus  the 
algorithm  runs  in  0(\V\  •  |f?[)  time.  For  a  placement  graph  it  is  not  difficult  to  show  that  both 
jFj  and  \E\  are  O(m),  and  thus  the  longest-paths  problem  can  be  solved  in  0(m2)  time  by  the 
general  algorithm.  But  a  linear-time  algorithm  can  be  found  by  exploiting  the  special  structure 
of  a  placement  graph  to  construct  a  list  £  of  length  O(m)  that  guarantees  the  correctness  of  the 
Bellman-Ford  algorithm.  We  now  look  at  the  structure  of  placement  graphs  more  closely. 

The  vertices  of  a  placement  graph  G  =  (V ,  E)  corresponding  to  the  chunks  on  the  top  of  the 
channel  have  a  natural  linear  order  imposed  by  the  left- to- right  order  of  the  chunks.  We  define 
the  partial  order  as  the  union  of  this  linear  ordei  with  the  similar  linear  order  of  bottom 

vertices.  Thus  u  v  for  vertices  u  and  v  if  their  chunks  lie  on  the  same  side  of  the  channel 

and  the  chunk  that  corresponds  to  u  lies  to  the  left  of  the  one  which  corresponds  to  v.  The 
left- boundary  vertex  Uq  precedes  all  other  vertices,  and  all  vertices  precede  the  right- boundary 
vertex  Vm  r  i-  The  partial  order  <  is  the  natural  extension  to  -<  that  includes  equality. 

'I  he  next  lemma  describes  some  of  the  structural  properties  of  placement  graphs.  Figure  10 
illustrates  the  impossible  situations  described  in  Properties  (i)  and  (it)  and  shows  the  only  kind 
of  simple  cycle  that  can  occur  in  a  placement  graph  together  with  the  two  consecutive  cross  edges 
that  satisfy  Property  (iii). 

Lemma  1.  Any  placement  graph  G  =  (V\/J)  has  the  following  properties: 

(i)  There  do  not  exist  cross  edges  (w,  v)  and  (x,y)  such  that  u  -<  x  and  y  -<  v. 

(u)  There  do  not  exist  cross  edges  (uf  u)  and  (x,y)  such  that  v  -<  x  and  y  -<  u. 

(in)  All  cycles  have  two  consecutive  cross  edges  (u,v)  and  (v,  tu)  such  that  w  u. 

f’roof.  Properties  (i)  and  (ii)  can  be  proved  by  considering  which  of  the  terminal  constraints 
from  Condition  (2)  induce  the  edges  in  the  placement  graph.  For  each  of  these  cases,  suppose  the 
edge  («,  t;)  was  caused  by  the  terminals  i  in  u  and  i  -f- 1  in  vt  and  the  edge  (x,  y)  came  from  the 


(a)  The  situation  forbidden  by  Property  (i). 


(b)  The  situation  forbidden  by  Property  («). 


(c)  Every  simple  cycle  contains  at  most  one  vertex  from  the  top  or 
at  most  one  vertex  from  the  bottom.  The  edges  incident  on  the 
vertex  are  a  consequence  of  Property  (in). 

Figure  10;  The  properties  of  the  placement  graph  enumerated  in  Lemma  1. 

terminals  j  in  x  and  j  +  t  in  y.  For  Property  (t)  we  have  u  -<  x  and  y  -<  u,  and  thus  i  <  j  and 
j *  <  i-M*  Canceling  t  from  this  latter  inequality  obtains  the  contradiction.  The  assumption 
to  be  proved  impossible  in  (it)  is  that  v  <  x  and  y  <  u>  which  implies  i  +  t  <  j  and  j  + 1  <  t. 
Since  t  is  nonnegative,  we  gain  a  contradiction. 

To  prove  Property  (in),  we  need  only  consider  simple  (vertex-disjoint)  cycles.  Since  no  cycle 
can  consist  solely  of  side  edges,  every  simple  cycle  must  have  a  cross  edge  (u,  v)  going  from  bottom 
to  top.  In  order  to  complete  the  cycle,  there  must  be  a  top- to- bottom  edge  (tu,  x)  such  that  v  <  w 
and  i  ;<  u.  If  v  =  w  or  x  =  it,  then  the  pair  of  edges  satisfies  Property  (iit).  But  if  v  ^  w  and 
then  the  pair  of  edges  violates  Property  (it).  | 

Each  edge  in  the  placement  graph  is  either  a  top  edge,  a  top- bottom  edge,  a  bottom- top 
edge,  or  a  bottom  edge.  For  each  of  these  four  sets  of  edges,  there  is  a  natural  linear  order  of 
edges  based  on  where  (u,  v)  precedes  (x,  y)  for  two  edges  in  the  same  set  if  u  <  x  and  v  <  y. 
Property  (ti)  guarantees  that  the  linear  order  holds  for  two  cross  edges  in  the  same  set.  Let  TT, 
TB,  BT,  and  BB  be  the  four  lists  of  edges  according  to  the  natural  linear  order,  and  include  the 
two  edges  out  of  t>0  and  the  two  edges  into  vm±\  in  cither  TT  or  BB  as  appropriate. 

The  list  C  used  by  the  Bellman-Ford  algorithm  is  constructed  by  a  merge  of  the  four  lists 
which  we  call  MERGE.  At  each  step  of  MERGE,  a  tournament  is  played  among  the  first  elements 
of  each  list.  If  (u,u)  and  (u,tu)  are  the  first  elements  of  two  lists,  then  (u,u)  beats  (v,u>)  if  tu  ^ 
u.  Since  there  may  be  more  than  one  edge  beaten  by  none  of  the  other  three,  ties  are  broken 


arbitrarily.  The  winner  is  appended  to  £  and  removed  from  the  head  of  its  JEt.  The  tournament 
i  i  then  repeated  until  no  edges  remain  in  any  of  the  four  lists.  The  performance  of  the  tournament 
can  be  improved  by  recognizing  that  only  six  of  the  twelve  possible  comparisons  of  edges  need 
be  tried,  and  that  w  ^  u  is  guaranteed  for  all  but  two.  Figure  11  shows  a  possible  ordering  of 
edges  in  £  for  the  placement  graph  in  Figure  8. 


In  order  for  MERGE  to  be  well-defined,  the  tournament  must  always  produce  a  winner, 
which  is  a  consequence  of  the  next  lemma. 

Lemma  2.  The  list  £  produced  by  MERGE  is  a  topological  sort  of  the  edges  of  E  according 
to  the  relation  R  where  (u,  v)R(vf  w)  tfw^u. 

Proof.  First,  we  show  that  the  relation  R  is  acyclic  so  that  the  edges  can  indeed  be 
topologically  sorted.  By  definition  of  R ,  a  cycle  in  R  induces  a  cycle  in  the  placement  graph. 
According  to  Property  (n't),  the  cycle  must  have  two  consecutive  cross  edges  (u,u)  and  such 

that  w  <  u.  But  since  (u,  v)R(v,  w),  we  also  have  that  w  ^  u,  which  is  a  contradiction. 

The  proof  that  MERGE  topologically  sorts  the  edges  of  E  according  to  R  makes  use  of 
the  fact  that  if  a  vertex  v  is  the  tail  of  an  arbitrary  edge  in  any  one  of  the  four  lists  TT,  TB, 
BT  or  BB,  then  for  every  u  <  v  there  is  an  edge  in  the  same  list  emanating  from  u.  Suppose 
that  MERGE  does  not  topologically  sort  the  edges  of  E  according  to  R.  Then  there  is  a  first 
edge  (u,  v)  in  £  such  that  there  exists  an  edge  (v,  it;)  earlier  in  £  and  (u,  v)R(vf  tu).  Consider  the 
edge  ( x ,  y)  in  the  same  list  as  (u,  t;)  that  competed  with  (vfw)  when  (v,w)  was  the  winner  of  the 
tournament.  For  each  of  the  possible  combinations  of  lists  for  (u,  v)  and  ( v,w ),  it  can  always  be 
deduced  that  there  is  an  edge  emanating  from  y  such  that  which  makes  (x,  y)  an  earlier  violator 
of  the  topological  sort  than  (u,  v).  | 

Since  each  edge  of  E  is  included  exactly  once  in  the  list  £  created  by  MERGE,  the  Bellman- 
Ford  algorithm  applied  to  £  has  a  running  time  linear  in  the  number  of  chunks.  The  correct  values 
for  longest  paths  are  produced  by  the  algorithm  if  for  every  vertex  u,  there  is  a  subsequence  of  £ 
that  realizes  a  longest  path  from  t/0  to  v,  under  the  assumption  that  there  are  no  positive- weight 
cycles  in  the  placement  graph.  Since  for  every  longest  path,  there  is  a  vertex-disjoint  longest 
path,  the  following  theorem  proves  the  correctness  of  this  linear- time  Bellman- Ford  algorithm. 

Theorem  3.  Let  G  be  a  placement  graph  with  left-boundary  vertex  vq.  Then  every  vertex - 
disjoint  path  beginning  with  Vo  is  a  subsequence  of  the  list  £  created  by  the  procedure  MERGE. 
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Proof.  We  need  only  show  that  eV’iry  pair  of  consecutive  edges  in  a  vertex-disjoint  path  from 
to  satisfies  R  because  then  Lemma  2  guarantees  that  the  path  is  a  subsequence  of  £ .  Suppose 
(•i,  t)  and  [v,  w)  are  two  consecutive  edges  on  a  vertex-disjoint  path  from  to  which  violate  R , 
that  is,  ti;  ^  u.  If  either  (u,t)  or  (t,tt)  is  a  side  edge,  the  pair  must  satisfy  R ,  and  thus  both 
must  be  cross  edges  with  the  vertices  u  and  w  on  the  same  side.  Since  if  w  =  t,  the  path  is  not 
vertex-disjoint,  we  need  only  show  that  w  -<  u  is  impossible. 

Assume,  therefore,  that  w  ■<  u>  and  consider  the  initial  portion  of  the  path  from  t’o  to  u . 
Since  t/o  -<  v  and  vo  ■<  w,  there  must  be  an  edge  (x,  y)  on  the  path  which  goes  from  the  set  of 
vertices  to  the  left  of  (i/,tt’)  to  the  right  of  (u,  u>)  in  order  to  get  to  u.  But  then  either  Property 
(t)  or  Property  (n)  is  violated  depending  on  whether  x  ^  v  and  w  -<  y,  or  y  -<  v  and  x  <  w.  | 


5.  Other  Wiring  Models 

The  reduction  from  the  fixed-separation  placement  problem  in  the  square-grid  model  to 
the  single-source-longest-paths  problem  is  possible  because  the  wirability  constraints  can  all  be 
written  in  the  form  vx  —  v3  >  6l3.  Thus  for  any  wiring  model  where  wiring  constraints  can  be 
written  in  this  form,  the  reduction  will  succeed.  Also,  it  should  be  observed  that  in  general,  the 
performance  of  the  single-source-longest-path  algorithm  will  not  be  linear,  but  will  be  a  function 
of  the  number  of  constraints  times  the  number  of  variables.  This  section  reviews  other  models 
and  gives  the  necessary  and  sufficient  wirability  constraints  for  each. 

1.  One-layer}  gridless  rectilinear  ([Dolev  et  al.}).  Wires  in  this  model  must  run  horizontally 
or  vertically,  and  although  they  need  not  run  on  grid  points,  no  two  wires  can  come  within  one 
unit  of  each  other.  The  wirability  constraints  for  this  model  are  the  same  as  for  the  square  grid 
model: 

d\ t  ^  t  and  ~ —  clx  ^  t 

f or  1  <  i  <  n  —  L  As  with  the  square-grid  model,  the  fixed-separation  placement  algorithm  for 
this  model  can  be  made  to  run  in  linear  time. 

2.  One-layer,  gridless,  rectilinear  and  forty-five  degree  ([Dolev  and  Siegel]).  This  model  is 
the  same  as  the  gridless  rectilinear,  but  in  addition  wires  can  run  which  have  slope  ^1.  The 
constraints  in  this  case  are 

t-r  —  b%  >  r\/ 2  —  t  and  bt+r  —  a%  >  r\pi  —  t 

for  tj\J 2  <  r  <  t  and  1  <  i  <  n  —  r.  The  placement  algorithm  for  this  model  runs  in 
0(min(tm2  -f  tn,  m3  +  In,  m3  -f  n2))  time. 

3.  One-layer,  gridless  ([Tompa]).  Wires  can  travel  any  direction.  The  constraints  are 

at+r  —  bx  >  y/r2 ~  t2  and  btq_f  —  at  >  Vr2  —  t2 
for  t  <  r  <  n  and  1  <  i  <  n  —  r.  The  placement  algorithm  runs  in  0(m3  +  n2)  time. 
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4.  Multilayer  models.  AH  the  modus  presented  until  now  have  been  one-layer  models.  It  is 
natural  to  generalize  to  Mayer  models  in  which  wires  may  travel  on  different  layers.  Remarkably, 
optimal  routability  can  always  be  achieved  with  no  contact  cuts  ((Baratz]),  that  is,  a  wire  need 
never  switch  layers.  The  necessary  and  sufficient  conditions  for  these  multilayer  models  are  a 
natural  extension  of  the  one-layer  conditions.  For  example,  in  the  one-layer,  gridless,  rectilinear 
model  the  conditions  are  modified  for  /  layers  to  be 

1-1  6,  ^  t  and  ——  u,  ^  t 


for  1  <  i  <  n  —  It. 

There  are  some  wiring  models,  however,  where  upper  and  lower  bounds  for  wirability  do 
not  meet.  For  these  models  a  constraint  graph  which  represents  upper  bounds  will  give  the  best 
possible  placement  for  those  bounds.  A  graph  representing  lower  bounds  will  give  lower  bounds 
on  the  best  possible  placement.  Together,  bounds  can  be  established  for  some  of  these  models, 
and  heuristic  algorithms  invoked  to  attempt  routing  within  the  feasible  range  of  optimality. 

6.  Extensions  and  Conclusions 

A  variety  of  related  placement  prob  ems  can  be  solved  by  the  method  described  in  this  paper. 
Some  entail  extensions  to  the  problem  specifications,  others  employ  different  wiring  models.  In 
this  section  we  shall  mention  a  few  extensions  we  can  handle  and  suggest  further  research  on 
more  complicated  problems. 

•  Nonriver  routing.  The  placement  algorithm  gives  optimal  placements  for  river  routing,  but 
there  are  other  routing  configurations  for  which  it  works  optimally  as  well.  One  example  is  the 
two-layer,  any-to-any  routing  problem  where  tv/o  sets  of  terminals  must  be  connected  across  a 
channel,  but  they  may  be  connected  in  any  order. 

•  Parallel  Channels.  Multiple,  parallel  horizontal  channels  are  easily  handled  within  the  same 
graph-theoretic  framework.  Each  row  of  chunks  is  represented  by  a  chain  of  vertices  (from  a 
common  left  margin  to  a  common  right  margin),  and  the  wiring  conditions  in  the  channels  are 
represented  by  edges  linking  adjacent  chains.  The  optimal  placement  is  achieved  by  solving  the 
whole  system,  using  the  Bellman-Ford  algorithm.  We  do  not  know  how  to  obtain  (or  show  the 
impossibility  of)  a  linear-time  algorithm  to  do  the  job  in  general,  in  contrast  to  the  special  case 
(one  channe1)  discussed  in  Section  4.  However,  the  number  of  edges  in  the  graph  is  still  linear  in 
the  number  of  vertices  (chunks),  thus  the  Bellman-Ford  algorithm  runs  in  time  0(n  +  m2)  (where 
n  is  the  number  of  nets,  m  —  the  number  of  chunks). 

•  Range- Terminals.  In  some  routing  situations  terminals  occupy  not  a  single  point,  but  rather 

a  contiguous  region  along  the  edge  of  the  channel.  For  example,  the  terminal  might  be  a  wire 
that  runs  along  the  edge  of  the  chunk,  and  connection  can  be  made  to  the  wire  anywhere. 
The  additional  flexibility  of  viewing  a  terminal  as  a  contiguous  range  of  points  can 

be  exploited  by  both  the  greedy  routing  algorithm  and  the  placement  algorithm  in  any  of  the 
river- routing  models  we  have  discussed. 

Each  range-terminal  is  specified  by  an  interval  [a^,a^J  or  [6^,6^].  The  greedy  routing 
algorithm  operates  as  before  with  minor  changes.  If  the  range-terminals  overlap,  the  wire  is 
routed  straight  across.  Otherwise,  assume  without  loss  of  generality  that  a **  <  and  use  the 
standard  greedy  algorithm  to  route  a  wire  from  a ^  to  6^. 
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The  wirability  conditions  for  placement  are  accordingly  adjusted.  In  th*  rectilinear  case,  for 
example,  the  condition  ax+t  —bx>tis  rewritten  as  a^_t  —  >  t  and  condition  bt+t  —  at>t 

becomes  b^t  —  a ^  >  t.  The  transformation  to  chunk  variables  is  as  before  and  the  placement 
algorithm  is  unchanged. 

An  interesting  extension  of  the  river-routing  problem  studied  here  is  the  two-dimensional 
problem  illustrated  in  Figure  12.  In  the  figure,  a  line  between  two  chunks  indicates  that  wires 
must  be  routed  between  them.  Unfortunately,  in  order  to  optimally  solve  this  general  problem, 
it  appears  that  the  constraints  indicated  by  the  lines  must  be  convex  in  both  dimensions,  not 
just  in  one  as  is  the  case  for  the  wiring  models  considered  here.  When  the  constraints  are 
convex,  however,  convex  programming  can  be  used  to  optimize  a  cost  function  such  as  the 
area  of  the  bounding  box  of  the  layout.  One  model  which  gives  convex  constraints  for  the 
general  two-dimensional  problem  is  the  one  in  which  all  wires  must  be  routed  as  straight  line 
segments  between  terminals  such  that  no  minimum  spacing  rules  are  violated.  This  model  is  not 
particularly  interesting  from  a  practical  standpoint,  however.  Heuristics  for  solving  the  related 
two-dimensional  compaction  problem  by  repeatedly  compacting  in  one  dimension  and  then  the 
other  can  be  found  in  (HsuehJ. 


Figure  12:  A  two-dimensional  extension  to  the  river-routing  problem.  A 
solid  line  between  two  modules  indicates  routing  occurs  between 
them. 

A  major  deficiency  of  placement  algorithms  is  that  they  lack  knowledge  about  the  wirability 
of  the  routing  problems  that  they  set  up.  We  have  shown  for  river  routing  that  wirability 
conditions  can  be  translated  directly  into  placement  constraints  without  the  overhead  of  wiring 
the  channel.  For  rectilinear  river  routing,  the  running  time  of  the  greedy  wiring  algorithm  is 
0(n2),  and  any  cost  function  for  placement  that  is  monotonic  in  spread  and  separation  can  be 
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o  .tLnized  in  0(n2)  time  without  the  overhead  of  routing.  Studying  wirability  in  the  general  case 
may  lead  to  the  development  of  heuristics  for  wirability  that  do  not  involve  routing.  A  program 
that  uses  this  heuristic  knowledge  should  be  able  to  outperform  the  iterative  place-route,  place- 
route  programs  that  dominate  today. 
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